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Abstract 


Measurements by Explorers 28, 33, 34 and 35 have been used to 
study the spatial characteristics of the magnetic field in the 
magnetosheath to a distance of 70 R E behind the earth. During 1966- 
1967 there were 1661 hours when one spacecraft of this multi-satellite 
system was monitoring the interplanetary medium while at least one 
other satellite was making magnetosheath measurements. This has 
made possible a separation of time and space variations and has 
permitted a study of the direction and magnitude of the magnetosheath 
field as a function of position and interplanetary field orientation. 
Results indicate that the magnetosheath field is several times the 
strength of the simultaneously measured interplanetary field in the 
sunward magnetosheath. This magnetosheath to interplanetary magnitude 
ratio decreases with distance from the subsolar point to values which 
are frequently less than unity at distances beyond 30 Rg and away from 
the bow shock. This ratio also displays a dawn -dusk asymmetry which 
is dependent on the interplanetary field orientation. T nterplanetary field 
lines perpendicular to the earth-sun line are associated with 
symmetrically distorted magnetosheach field lines in the dawn and 
dusk hemispheres and are consistent with the draping of field lines 
around the magnetosheath „ When the interplanetary field is aligned 
near the spiral angle, fields measured in the dusk hemisphere are 
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much more ordered than those measured in the dawn hemisphere behind 
the earth. These experimental results are in general agreement with 
the theoretical predictions of the magne togas dynamic models. 
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Introduction 

4 

Satellite measurements indicate that solar magnetic fields frozen 
in the solar plasma are connected away from the sun by the solar wind 
(Ness and Wilcox, 1966; Fairfield, 1968). The problem of solar wind 
plasma flow past the geomagnetosphere has been studied theoretically by 
Spreiter et al. (1966,1968), Alksne (1967), and. Dryer and Heckman (1967) 
using a magnetogasdynamic approach. The plasma flow is determined 
from a gasdynamic solution which is then used to calculate the three 
dimensional deformation of the geomagnetic field in the magnetosheath. 
The predicted spatial features of the magnetic field may be summarized 
as follows: 

1. The ratio of the magnitude of the magnetosheath field to 
the simultaneous interplanetary field, B/B q , varies from 
values greater than 4 near the stagnation point to values 
less than unity in the downstream magnetosheath; 

2. Interplanetary field lines convected into the magnetosheath 
become draped about the magnetosphere in a manner 

that depends on the direction of B q . This dependence results 
in a dawn-dusk asymmetry of both direction and B/B 0 
distribution for oblique orientations of the external field. 

Experimental measurements with the IMP 2 satellite in the noon- 
dawn region of the magnetosheath showed a general agreement with 
these predictions (Fairfield, 1967). The results indicated that the 
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magnetic field is convected through the earth's bow shock where it 
undergoes an increase in magnitude but relatively little change in 
orientation. Fairfield further concluded that as the field is 
convected deeper into the magneto sheath, it is subjected to a greater 
distortion in direction until it is aligned tangent to the magnetopause. 
Now a considerable number of additional measurements are available to 
extend the experimental study to the region of the magnetosheath behind 
the earth on both the dawn and the dusk sides of the magnetosphere. 

In practice the steady-flow constant field conditions assumed 
in the theory do not exist and the normal time variations of the 
interplanetary field are also seen in the magnetosheath. These time 
variations present the primary difficulty in comparing experimental 
results and the time independent theory. The best technique for 
overcoming this difficulty is to use simultaneous satellite measurements. 
Since the typical scale of interplanetary field variations is large 
compared to the dimensions of the magnetosphere, it is reasonable to 
consider an interplanetary measurement of B 0 on any earth-orbiting 
spacecraft as representative of conditions across the entire extent 
of the bow shock. It is then possible to investigate simultaneous 
data from a second spacecraft measuring 5? in the magnetosheath and to 
attribute the difference in vector measurements to interaction with 
the shock and magnetosphere. This technique has been utilized in 
analyzing 1661 hours of simultaneous data obtained by the satellites 


* 6 “ 

Explorers 28, 33, 34 and 35 during the period July 1966-January 1968. 
The spatial distribution of B/B 0 and the directional characteristics 
of the magnetosheaf> field have been determined in the region back 
to 70 behind the earth. These results have been compared to the 
theoretical results and general agreement has been found. 
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Satellites and Instrumentation 

The Explorer 28 (IMP 3) satellite was launched on May 29, 1965 
Into a highly eccentric earth orbit with an initial apogee of 
approximately 42 Rj? and an orbital period of 5.8 days. A single mono- 
axial magnetometer sampled the vector magnetic field once every 40.96 
seconds with digitization error of £0.4 y . The instrumentation and 
data analysis were similar to those of IMP 1 (Ness et al . , 1964) and 
IMP 2 (Fairfield and Ness, 1967). 

The Explorer 33 satellite was launched into an extremely high- 
apogee high -perigee orbit on July 1, 1966. This spacecraft had an 
initial apogee of 70 Rg and orbital period of 14 days. A triaxial 
fluxgate magnetometer, with a sensor range of +64 y, measures the 
vector magnetic field every 5.11 seconds with a digitization error 
of +0.25 y . The instrumentation has been described in detail by 

f a 

Behannon (1967). 

Explorer 34 was placed in a high inclination (67.4°) orbit on 
May 24, 1967, with an apogee of 34 Rg and an orbital period of 4.3 
days. A. triaxial fluxgate magnetometer makes a vector magnetic 
field measurement every 2.56 seconds with a digitization error of 
+0.16 Y and +0.64 Y in the low and high sensor ranges of +32 Y and 
+128 Y , respectively. The instrumentation has been described by 
Fairfield (1968). 
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The Explorer 35 satellite was launched on July 19, 1967 and 
injected into lunar orbit on July 22, This spacecraft performs 
magnetic field measurements of the Interplanetary field, the 
magnetosheath field and the geomagnetic tall field at an approximately 
constant geocentric distance of 60 R^. The experiment and analysis 
are essentially identical to those used on Explorer 33 except that 
a range switch was added to the Explorer 35 instrument which 
automatically choses the appropriate £24 7 or £64 V sensor range with 
its corresponding £.09 7 or .25 7 digitization errors (Ness et al., 
1967). 
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Analysis 

Figure 1 shows the relative positions of Explorers 33, 34 and 35 
in September -October 1967. On September 20-21, 1967, Explorer 34 was 
at apogee in the interplanetary medium on the sunward side of the earth 
while Explorer 33 was in the dusk magnetosheath and Explorer 35 was 
in the dawn magnetosheath. Hie relative positions of these spacecraft 
at other times during the interval of interest has been published 
elsewhere (Behannon et al. 1968). 

Figure 2 shows typical detailed simultaneous magnetosheath and 
interplanetary data (81.81 second averages for Explorer 33 and 20.45 
second averages for Explorer 34) for a 12 hour period on September 20- 
21, 1967 when the two spacecraft were separated by , distance of 
approximately 75 R E along the earth -sun direction. Figure 2 shows 
that the magnetic field directions at the two satellites, as given 
by the solar ecliptic latitude angle 6 and azimuthal angled , were 
quite similar. The field magnitude trends were also similar; however, 
the magnetosheath field magnitude was enhanced relative to the 
interplanetary magnitude. 

Large field discontinuities such as thos# near 2200 UT arrive 
at the two spacecraft with time delays which are consistent with 
fr< zen-in magnetic fields being convected with the plasma velocity. 

Even for this case of near maximum separation of the two spacecraft, 
the time delay between observations of field discontinuities at the 
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respective satellites is only of the order of 15 minutes. Since this 
time delay is always small in comparison to one hour and since the 
field usually changes by only a small amount from one hour to the 
next, averages for the same hour have been considered as representative 
of "simultaneous " interplanetary and magnetosheath conditions. Hourly 
average magnitudes were computed from the hourly averages of the 
field components. 

Figure 3 shows typical hourly average field magnitudes and 
directions, in this case for Explorer 33 in the dawn magnetosheath 
and IMP 3 in the interplanetary medium on the sunward side of the 
magnetosphere. Five days of simultaneous observations by the two 
spacecraft from an outbound passage of IMP 3 through the bow shock 
(0500 UT August 27) to the next inbound passage through the shuck 
(2000 UT September 1) are shown. Field magnitudes and directions 
generally are not identical but good correlation between changes in 
the fields is seen. The observed gradual increase with time in the 
magnetosheath field -magnitude enhancement is due to the movement of 
Explorer 33 toward the bow region of the magnetosheath. The occurrence 
of transient disturbances such as those associated with the magnetic 
storm sudden commencements at 1315 UT on August 29 and 1112 UT on 
August 30 emphasizes the need to normalize the magnetosheath field 
magnitude to the simultaneous interplanetary field magnitude in order 
to separate the spatial variations. 


- 11 - 


To facilitate comparison of experiment and theory, both the 
satellite positions in the magnetosheath and the field vectors have 
been transformed from solar ecliptic coordinates to the coordinate 
system used by Spreiter et al. (1968). This transformation consists 
of a rotation of less than 90° about the solar ecliptic X axis until 
the Y axis is in the plane defined by the X axis and the interplanetary 
magnetic field vector. In this plane interplanetary fields oriented 
in the spiral angle quadrants 90° < 9se < 180° and 270 o<c Pse < 360° 
make angles 0 o< <90° with the X axis. Fields oriented in the 

quadrants 0 o< cp gE < 90 or 180° < 9 gE < 270° make angles with the X 
axis of -90 o< t|i <0°. This is illustrated in Figure 4. In this new 
coordinate system the geomagnetic dipole will assume a variety of 
orientations but if the magnetosphere is a blunt body exhibiting 
symmetry about the X axis, as is assumed in the theories of hydromagnetic 
magnetospheric flow, this will not be an important consideration 
and the new coordinate system will be appropriate. 

The 1661 hourly average magnetosheath field -magnitude values 
were sorted according to the corresponding interplanetary field 
orientations using eighteen 10° \jf angle sectors. Figure 4 shows the 
resulting data distribution as a function of \|r . The asymmetry in 
the distribution reflects both the preference of the interplanetary 
field for the spiral angle and a three dimensional geometrical effect. 

The latter effect is due to the fact that the earth-sun line is a 
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pole about which all possible orientations of the plane of the 
interplanetary field form a spherical distribution of field vectors. 

In this distribution the number of possible orientatiors of the 
field increases away from the pole. 

The selection of magnetosheath data for this study was carried 
out with particular care. Generally hours during which a spacecraft 
was sampling only magnetosheath fields could be identified unambiguously. 
On those occasions when traversal of either the magnetopause or the 
bow shock wave was not clear, data was omitted rather than risk 
inclusion of magnetosphere or interplanetary data. 
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Magnitude Distribution 

The general agreement of the observations with theoretical 
contours of B/B 0 is illustrated in Figure 5. The dashed theoretical 
contours shown are reproduced from Alksne (1967) and were computed 
for a plane offset from the XY plane by approximately 6.7 R^, using 
a Mach number M=8 and a specific heat ratio of 5/3. Theoretical 
contours are labeled with circled numbers and experimental measurements 
of B/B 0 are indicated by uncircled numbers. The measurements were 
performed at distances from the XY plane lying within the | Z | range 
of 4-15 Rg. Some values have been omitted because of crowding. 
Theoretical data correspond to an interplanetary magnetic field 
orientation of V =90° while experimental data are for interplanetary 
fields with angles 80° < V < 90°. For this case a magnitude 
distribution pattern which is symmetric with respect to the earth-sun 
line is predicted. 

In all cases studied, the experimental data are in general 
agreement with the theory even though the data used correspond to a 
wide range of Z coordinates^ a 10° range of ifi angles, and varying flow 
parameters for the solar wind during the long period covered by the 
observations. The measured magnitude ratios in Figure 5 are seen to 
decrease from values as high as 6 near the bow to values less than one 
downstream. A predominance of values less than one are found at 
greater distances downstream. 
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Although 1661 hourly averages were used in this study, they 
were spread over a large region of space and were divided among 
18 angular intervals in the analysis. Eighty-seven percent of these 
hours are from positions in the magnetosheath with X<-10 Rg, a region 
not covered in detail by the published theory. For these 
reasons, experimental contours could not be drawn with the resolution 
of the theoretical contours. In spite of these limitations, statistical 
analysis of the data yields results in agreement with the principal 
features predicted by the theory. 

Figure 6 shows averages of the magnitude ratio computed for 
intervals of 10 along the earth-sun line from 20 Rg sunward of 
the earth to 80 R E in the antisolar direction. Averaged together 
are magnetosheath data from all values of Y and Z in each X interval, 
and for all orientations of the external field. The error bars 
represent the standard deviations computed for each 
10 R„ interval. The continuous decrease is seen from a gross average 

a 

of approximately 4 to values of the order of one and finally to less 
than one well downstream. There are many cases of B/B 0 < 1 beyond 

-50 R but higher values near the shock tend to raise these averages. 

E 

The gradient can be described as an exponential decrease as is shown 
in Figure 6. 

More complex asymmetrical effects can be seen by dividing the data 
into dawn and dusk distributions for various interplanetary field 
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orientations. There were 1082 cases of the interplanetary direction 
lying in the dawn quadrant containing the spiral angle, and 579 cases 
in the dusk quadrant. The results of this are summarized in Figure 7. 
These data demonstrate, in a statistical sense, the dependence of 
magneto sheath magnitude on the orientation of the external field. 
Distributions for the interplanetary field orientation in the quadrant 
of the spiral angle (A) and in the perpendicular quadrant (B) are 
shown. Theory predicts that for orientations in the spiral angle 
quadrant the magnetosheath field should be stronger on the dusk (+Y) 
side than on the dawn (-Y) side. This relationship is generally 
seen for these distributions. For orientations in the dusk quadrant 
the situation should be reversed, and except for one interval that 
was found to be the case. 
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Magnetoshcath Field Directions 

To study the spatial variations of magnetic field direction in 
the magnetosheath the hourly average field vectors were normalized 
in magnitude to the simultaneous interplanetary field magnitudes and 
projected on the XY plane. Two cases, corresponding to 80° < 90° 

and 30 < f <40 , arc shown in Figures 8 and 9, respectively. For 
simplification of presentation, field polarities in the magnetosheath 
were reversed whenever the corresponding interplanetary measurements 
had a component directed from the dawn to the dusk meridian. The field 
vectors in Figures 8 and 9 that are located near the X axis correspond 
to measurements made at large jz |. The direction of these vectors 
is generally found to be near the interplanetary field direction. 

One of the most striking features observed in the data for the 
various interplanetary field orientations was a marked lack of 
symmetry in the ordering of the vectors on the dawn and dusk sides. 

When the field is aligned nearly perpendicular to the earth-sun 
line as in Figure 8, the vectors are quite orderly and consistent 
with the draping pattern back to X » - 70 R E . For other interplanetary 
field orientations where there is a smaller angle between the field 
and the earth-sun line, a smaller percentage of the orientations are 
consistent with this draping pattern in the dawn hemisphere. Figure 
9 is typical in that a rather wide variety of orientations are seen 
in the dawn hemisphere while the dusk hemisphere remains quite orderly. 
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This difference between the hemispheres may be attributed to 
conditions ah the shock when it is realized that the interplanetary 
field is oriented nearly parallel to the average bow shock position 

in the dusk hemisphere but is approximately perpendicular to the shock 

<1 

in the dawn hemisphere u Dusk side fields are then able to convert 
through the shock with little angle change but dawn side fields must 
undergo a large angle change in order to conserve the normal component 
of magnetic field flux across the shock. When the field is oriented nearer 
the earth-sun line, the minimum angle change necessary to conserve 
flux is achieved by rotating the field in the direction opposite to 
that expected from draping . This effect can be seen in the theoretical 
work (Alksne, 1967). 

To determine whether the magnetosheath field lines were bent clock- 
wise (looking down from the north pole) or counterclockwise relative to 
the interplanetary orientation, vector plots in the dawn hemisphere 
were studied. The number of cases clearly in agreement with the draping 
orientation and those in agreement with the opposite tendency were 
tabulated and are shown in Table 1. The last column contains the number 
of magnetosheath measurements in the dawn hemisphere associated with 
the ranges of interplanetary ij; values listed in. the first column. The 
intermediate columns 2 through 5 exhibit the separation of these cases 
according to whether the magnetosheath measurement is rotated clockwise 


or counterclockwise. 
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Columns 6 and 7 list cases where the angular difference between 
the interplanetary field and the magnetosheath field projection in the 
XY plane was less than 10°. Only the XY field components were considered 
in preparing Table 1, even though Z components were sometimes large. 

Large Z components tended to be associated with the cases tabulated 
in columns 4 and 5. 

Results indicate that fields perpendicular to the earth-sun line 
are almost always rotated in the direction consistent with both the 
draping of field lines and minimum deviation of the field at the shock. 
For smaller i|r angles the percentage of fields rotated in this direction 
decreases. For tjf <40° where the conditions of draping and minimum 
angle change dictate opposite directions of rotation, there are an 
approximately equal number of rotations in either direction. Because 
most of the experimental data used correspond to large negative X 
positions, the field direction results cannot be directly compared 
with the theory since detailed theoretical results are not available 
for that region of the magnetosheath. 
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Summary of Results 

Simultaneous observations from four satellites have produced 
an experimental mapping of the magnetic fields in the magnetosheath 
that is in zeroth order agreement with the predictions of magneto- 
gasdynamic theory* These experimental results support the theoretical 
prediction that the ratio of the magnetosheath field magnitude 
relative to the interplanetary field magnitude, has a value of 4 or more 
near the bow but decreases to less than unity downstream. The 
observations also support the predicted dependence of the magnetosheath 
field magnitude at a given point on the orientation of the field 
external to the magnetosheath* In particular the analysis reveals a 
dawn-dusk asymmetry in the magnitude distribution in which the 
magnitude difference between dusk and dawn reverses sign as the inter- 
planetary field orientation changes from the spiral angle direction 
to the spiral angle + 90°. 

The experimental results also support the predicted dawn-dusk 
asymmetry in field orientation* When the interplanetary field is 
aligned near the spiral angle, the hourly average field vectors in 
the dusk hemisphere are more ordered than those in the dawn magneto- 
sheath* Distortion of the field direction due to draping extends at 
least as far downstream as -70R E in the case of the external field 
perpendicular to the earth-cun line* The tendency for the downstream 
field to have an orientation consistent with draping decreases as 
the angle between the interplanetary field and the earth-sun line 
decreases, suggesting that conditions at the bow shock may sometimes 
be the more important factor influencing the magnetic field direction. 
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Although the experimental data presented above clearly agrees 
with the steady state magnetogasdynamlc theory in the zeroth order, 
there is still a considerable amount of scatter among the hourly 
average data points. Some of this scatter is undoubtedly due to 
the inadequacies in removing time variations from the data. Errors 
are involved in ignoring convection time delays when comparing hourly 
averages and in assuming interplanetary measurements as representative 
of fields across the entire shock surface. Scatter also arises from 
comparing measurements made at different positions and at times when 
the values of various solar wird parameters are different. In 
addition to these experimental difficulties, it is possible that 
first order effects included in the steady-state magnetogasdynamic 
theory are responsible for some scatter. Geomagnetic dipole 
orientation might be important in a first order treatment and the 
assumption of axial symmetry used in the zeroth order theory would 
no longer be valid. 
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TABLE 1 
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Figure Captions 


Figure 1 


Figure 2 


Figure 3 


Figure 4 


Figure 5 


Solar ecliptic XY plane projections of trajectories of 
Explorers 33, 34 and 35 during September 17 -October 16, 
1967. Orbits 27-31 are shown for Explorer 34. Distances 
are in earth radii (Behannon, Fairfield and Ness, 1968). 
Simultaneous magneto sheath and interplanetary magnetic 
field measurements from Explorers 33 and 34. The field 
magnitude F and solar ecliptic latitude and longitude 
angles 6 and cp are shown. 

Simultaneous hourly average magnetic field data from 
Explorer 33 in the dawn side magnetosheuth and IMP 3 in 
the interplanetary medium during August 28-September 1, 
1966. 

Distribution of magnetosheath data used in the analysis 
as a function of interplanetary field orientation. As 
is shown schematically, t is the angle between the 
interplanetary field and the earth-sun line in the plane 
of the field vector and the earth-sun line. 

Comparison between experimental values of field magnj tude 
radio B/B q and theoretical contours, for the case of the 
interplanetary field oriented transverse to the earth-sun 


line. 
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Figure 6 


Figure 7 


Figure 8 


Figure 9 


Averages of normalized magneto sheath magnetic field 
magnitude for intervals of 10 R-g along the earth-sun 
line showing the continuous decrease of the magnetosheath 
field with distance downstream to values less than the 
interplanetary field magnitude. 

Illustration of asymmetry in average field magnitude 
between dawn and dusk sides of the magnetosheath. When 
the interplanetary field is aligned in the quadrant of 
the spiral angle (case A) the field magnitude on the dusk 
side is shown generally to be larger than for orientations 
corresponding to case B. The reverse is seen to be true 
for the dawn side field. 

Projection of normalized magnetosheath field vectors onto 
the XY plane defined by the interplanetary field and the 
earth-sun line, for cases when the interplanetary field 
was approximately transverse to the earth-sun line. Note 
the high degree of symmetry between the field patterns 
on the daTfn and dusk sides. 

XY plane projection of magneto sheath field vectors for 
cases when the interplanetary field was oriented at 30° - 
40° from the earth-sun line. In this case less symmetry 
between dawn and dusk fields is evident than is seen in 


Figure 8. 
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